METHOD AND APPARATUS FOR BALANCING 



This application is a continuation of copending nonprovisional application 
5 number 09/282,755 filed on March 31, 1999. 

BACKGROUND OF THE INVENTION 
(1) FIELD OF THE INVENTION 

10 

This invention generally relates to a method and an apparatus for balancing and 
more particularly to a method and to an apparatus for balancing a rotating machine tool 
and/or rotating/moving machine tool assembly. 

15 (2) BACKGROUND OF THE INVENTION 

Machine tool assemblies are used to create and/or to form various parts, products, 
and other types of entities and/or apparatuses. Particularly, these apparatuses include 
tools which are usually and removably placed, by means of a tool holder, upon a rotating 
or moving spindle and which are adapted to engage a workpiece or other apparatus and to 

20 engagingly form the workpiece into a desired shape and size. While these machine tool 
assemblies adequately create the desired parts, products, and/or other entities they are 
prone to imbalance failure and/or imbalance error. Hence, the precision and reliability of 
these relatively high-speed machine tool assemblies is limited by the imbalance which 
exists and/or arises as the assembly is used (e.g. as the tool is engaging rotated and/or 

25 moved against the workpiece). Similar types of unbalance or imbalance difficulties arise 
in the rotation and the utilization of relatively high-speed turbomachinery which is used 



in a variety of manufacturing and electrical power generation facilities and which causes 
the loss of significant revenues if "shut down" and/or disabled for even a relatively brief 
period of time. 

Some balancing devices have been created and utilized in order to overcome these 
5 imbalance difficulties. While somewhat effective, these prior balancing devices and 
methodologies suffered from various drawbacks. For example and without limitation, 
many of these prior balancers required that the tool assembly and/or rotating machinery 
be stopped or "interrupted" before a balance was achieved. This "interruption" or "off 
line" balance was and remains relatively costly and highly inefficient. Some of the other 

10 prior balancing devices and methodologies experienced and/or caused a relatively high, 
although oftentimes temporary, increase in the overall vibration or imbalance level as the 
balance "correction" was implemented. This transitory high vibration level was usually 
caused by the movement of the balancing portion of the assembly to a new balancing 
position and oftentimes caused the tooling assembly to be imprecise and to fail even 

15 before overall balance was achieved. The very act of balancing the assembly causes the 
assembly to fail or perform in an undesirable manner. Further, many of these prior 
balancers did not and do not adjust the imbalancing "thresholds" (e.g. the vibration 
values at which balancing was and/or is to begin and end) during the "working life" and 
for various speeds of the machine. Hence, as the machine (tooling assembly) ages and as 

20 the machine and/or tooling assembly employs a variety of speeds its needs for balancing 
changes. Unfortunately, these new and different balancing needs were not and are not 
properly reflected by these prior deployed balancer assemblies. 
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There is therefore a need for a method and a balancing apparatus to allow for the 
relatively efficient balancing of a machine tool assembly, a relatively high-speed 
turbomachine assembly, or virtually any other rotating and/or moving machine and/or 
apparatus and which substantially prevents temporary vibration level increases during 
5 balancing and which allows balancing thresholds to be dynamically modified over the 
"working life" of the machine. Applicant's invention addresses these needs and these 
aforedescribed drawbacks of these prior art balancing devices and methodologies in a 
new and novel fashion. Importantly, Applicant's invention is not limited to the balancing 
of a machine tool assembly. Rather, Applicant's invention is directed to the balance of 
1 0 any moving and/or rotating member. 

SUMMARY OF THE INVENTION 

It is a first object of the invention to provide a method and an apparatus for 
15 balancing. 

It is a second object of the invention to provide a method and an apparatus for 
balancing which overcomes some or all of the various drawbacks of the prior balancing 
apparatuses and methodologies, including for example and without limitation, those 
which have been previously delineated above. 
20 It is a third object of the invention to provide a method and an apparatus for 

balancing which overcomes the various drawbacks associated with prior balancing 
apparatuses and methodologies and which allows for the dynamic balancing of a tool 
assembly and/or other types of rotating machinery and apparatuses. 
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According to a first aspect of the present invention a balancer is provided. The 
balancer, according to this first aspect of the present invention, includes a balancer to 
selectively balance unbalance in a rotating tool assembly; and control means, connected 
to the balancer, for causing said balancer to calculate an influence coefficient of the 
5 rotating tool assembly and to thereafter selectively balance the rotating tool assembly by 
the movement of the balancer, movement being defined by said calculated influence 
coefficient. 

According to a second aspect of the present invention a balancer is provided. The 
balancer includes a balancing means to selectively balance a rotating tool assembly; and 

10 control means, coupled to said balancing means, for dynamically calculating vibration 
levels necessary for said balancing to begin and to end. 

Further objects, features, and advantages of the present invention will become 
apparent from a consideration of the following detailed description of the preferred 
embodiment of the invention, the claims, and the attached drawings. 

15 According to a third aspect of the invention, a balancer methodology is provided. 

This methodology includes the steps of a method of balancing a tool assembly by 
selectively employing a balance weight correction to said tool assembly, said method 
comprising the steps of employing a first balance weight correction to said tool assembly, 
measuring an amount of vibration associated with said balance tool assembly, estimating 

20 an influence coefficient value of said tool assembly, dividing said measured amount of 
vibration by said estimate influence coefficient thereby creating a certain value; 
multiplying said certain value by a gain parameter value, thereby creating a second 
certain value and subtracting said second certain value from said first balance weight 
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correction, thereby creating a new balance weight correction and applying said new 
balance weight connection to said tool assembly. 

According to a fourth aspect of the present invention, a balancer is provided 
which substantially reduces and/or eliminates transitory vibration increases. 
5 Further objects, features and advantages of the present invention will become 

apparent from a consideration of the following detailed description of the preferred 
embodiment of the invention, these included and subjoined claims and the attached 
drawing. 

BRIEF DESCRIPTION OF THE DRAWINGS 

10 The following drawing is useful in understanding the nature and the operation of 

the preferred embodiment of the invention. The drawing is as follows: 

Figure 1 is a block diagram of a balancing apparatus made and operating in 
accordance with the preferred embodiment of the invention and shown operatively 
coupled to and deployed upon a typical tool assembly. 

15 

DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENT 

OF THE INVENTION 

Referring now to Figure 1 there is shown a balancing assembly 10 shown in 

20 operative and assembled relation with a typical .and conventional tool assembly 12. 

Assembly 10, as shown, is made in accordance and operates in accordance with the 

preferred embodiment of the invention. 

As shown, tool assembly 12 includes a moving or selectively rotatable spindle 14 

having an integral tool holder 16 into which a tool 18 is removably placed. In operation, 

25 the movement and/or rotation of spindle 14 causes the tool 18 to move and/or rotate and 
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to engage a workpiece 20 for the purpose of "machining" the workpiece 20 into some 
sort of desirable shape, size, and/or geometry. As earlier explained, tool assembly 12 is 
subject to unbalances which cause the tool 18 to perform imprecisely and undesirably and 
which causes the tool 18 and/or other portions of assembly 12 to fail and/or to become 
5 fatigued. 

As shown, balancer assembly 10 includes a controller 22 which in one 
embodiment of the invention comprises a microprocessor acting and/or operating under 
stored program control and an electrical driver 24 which is selectively coupled to the 
source of electrical power 26 through controller 22 and which is adapted to selectively 

10 couple electrical power to the balancer actuators and/or rotors 28, 30 in a manner 
effective to correct a measured and/or calculated unbalance condition associated with the 
tool assembly 12, by movement of these rotors 28, 30. Such unbalance measurement 
may be made by the use of one or more vibration sensors 32. Moreover, balancer 
assembly 10 further includes a position sensor 34 which is communicatively and 

15 operatively coupled to the controller 22 and which is adapted to provide positioning 
information associated with the placement and/or positioning of the rotors 28, 30, to the 
controller 22. Applicant's invention is directed to the manner of moving these unbalance 
correction rotors or other balance correction members. Thus Applicant's invention may 
be used with a wide variety of computerized balancing systems; each of these systems 

20 being adapted to selectively provide a balancing weight correction achieved by the 
movement of a certain member or members each having a certain non- symmetrical or 
uneven mass distribution. One example of such a balancing system 10 is shown for 
example and without limitation, within United States Patent Number 5,757^662; which 
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was filed on November 29, 1994; which issued on May 26, 1998; which is assigned to 
Applicant's assignee; and which is fully and completely incorporated herein by reference, 
word for word and paragraph for paragraph. Another non-limiting example of a 
balancing system 10 is shown within United States Patent Number 4,977,510; which was 
5 filed on July 21, 1989; which issued on December 11, 1990; which is assigned to 
Applicant's assignee; and which is fully and completely incorporated herein by reference, 
word for word and paragraph for paragraph. Yet another balancing and control system is 
shown within United States Patent Number 5,231,265; which was filed on September 28, 
1990; which issued on July 27, 1993; which is assigned to Applicant' s assignee; and 

10 which is fully and completely incorporated herein by reference, word for word and 
paragraph for paragraph. The use of Applicant's improved, new, and novel balancer 
control methodology will now be explained in further detail. 

According to Applicant's preferred balancing methodology, as shown below, an 
automatic and/or substantially automatic tuning mechanism is provided which is based 

15 upon two adaptive parameters which are calculated during each control iteration based 
upon the convergence of the influence coefficient estimate. These parameters are used to 
define the correction movement of the balancer 10 and/or of those movable balancer 
members of the balancer assemblies in order to balance a tool assembly spindle, 
turbomachinery, or other moving member. 

20 The measured vibration of the tool assembly 12 can be considered to be an "error 

signal" that is desired to be eliminated by the use of an unbalance compensation 
apparatus, such as and without limitation balancer assembly 10. In the preferred 
embodiment of the invention, during each control iteration, vibrational data is sampled 
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and depending upon the "level" or amount of reported unbalance, the rotors 28, 30 are 
moved to a certain position. That is, vibrational data emanating from the sensor(s) 32 is 
communicatively coupled to the balancer controller 22 where it is "recognized" and 
where it causes certain balance operations to occur in order to minimize the amount of 
5 unbalance present within the tool assembly 12. This "error signal" can be mathematically 
represented as a complex phasor having a certain phase angle (e.g. obtained by the use of 
a fixed target on rotating assembly which is utilized by the position sensor 34 ) and a 
magnitude. 

The total vibrational error signal can be thought of as the linear and mathematical 
10 summation of all of the system "disturbances" (e.g. anything of a physical or electrical 
nature that might cause the tool system 12 to become unbalanced, such as and without 
limitation a "wearing away" of a portion of the tool holder) and of the overall effect of 
the balance corrections which were made to the system. Hence, Applicant has discovered 
that the total system vibrational error "e k " can be computed as follows: 
15 e* = C(co)W k +d Eq.(3.1) 

Where "C(a>)" is the influence coefficient of the tool assembly 12 and 
Where "d" is the total disturbances applied to and/or resident within tool assembly 12 and 
where "Wk" is the balance weight correction applied to assembly 12 by the balancing 
assembly 10. 

20 Hence, according to the relationship set forth in Equation 3.1^ the vibrational error 

associated with and/or "reported'Vcalculated/measured at of each control iteration is the 
mathematical sum of all of the disturbances present during or at that particular iteration 
and the mathematical multiplicative product of the balance correction made during that 
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interval with the influence coefficient of the machine. In essence, as should now be 
apparent to those of ordinary skill in the balancing art, the influence coefficient "C(co)" is 
a measure of the total response of assembly 12 to a certain amount of unbalance 
correction and is calculated and/or estimated in measurement units of units of vibration 
5 per unit of unbalance. Essentially, the metric "C(co)" defines and/or is a measure of the 
response of the system 12 to a certain amount of applied unbalance compensation and is a 
transfer function. Moreover, as the assembly ages the value of this influence coefficient 
may change. 

Assuming that the disturbance remains substantially similar across many control 
10 iterations and assuming that the amount of vibrations in the "k+1" period is relatively 
small, Applicant has discovered that one can develop a mathematical relationship 
between corrections in succeeding correction intervals. This is shown below: 

W k+1 =W k -^/c Eq. (3.2) 

15 where c is an estimate of the actual influence coefficient c. 

Applicant has discovered that the use of a gain parameter "oT enhances stability 
and when used in combination with Applicant's preferred embodiment of the invention 
reduces and/or substantially eliminates transitory balancing vibration increases. This is 
shown below: 

20 

e k 
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Eq. (3.3) 

This gain parameter a can improve control robustness to errors in the influence 
estimate c. We can, without loss of generality, consider the disturbance d to be 
represented by the influence of a hypothetical unbalance u (which may include any 
5 synchronous forcing, not just unbalance) acting through some cumulative influence 
coefficient c u such that 

d=cc u u Eq. (3.4) 

Equation (3.1) can be written 

ek=c(w k +c u u) Eq. (3.5) 
10 We can define a cumulative forcing term where 

v k ^w k +c u u Eq. (3.6) 

Substituting Eq. (3.6) into Eq. (3.5) gives 

e k =cv k Eq. (3.7) 

Again assuming that the disturbance u and influence coefficients c and c u do not change 
15 over one control iteration, we can use Eq. (3.7) and rearrange Eq. (3.6) for both control 
iterations k and k+\ in Eq. (3.3) giving relationship: 

n + i = n-a7V, Eq.(3.8) 
c 



v M - = (l-«f)v A Eq.(3.9) 



The recursive forcing of Eq. (3.9), combined with Eq. (3.7) yields the controlled vibration 
20 error history 



C 



e k+ t= 1-^7 e o Eq.(3.10) 
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It then follows that the single-plane active balancing control is stable for constant 
influence coefficients c if and only if 



R = 



c 
c 



<1 Eq.(3.11) 



5 This criterion can be represented graphically in the complex plane as shown in Figure 
3.2. 

From Figure 3.2, one can see how the control gain parameter a can enhance 
robustness. For the case in which the actual influence coefficient does not change during 
control convergence, it can be shown, based on Eq. (3.11), that the non-adaptive control 
1 0 law of Eq. (3.3) is stable if and only if 

a<2^co^6 c -e £ ) Eq. (3.12) 

where \c\ and \c\ are the magnitudes of the actual and the estimated influence coefficients 
and 0 C and 0 e are their respective phase angles. Since a is defined to be greater than 

zero, if the phase angle of the influence coefficient estimate lies between ±90° 
15 (noninclusive) of the actual influence coefficient phase angle (i.e., the real part of the 
quotient c/c is greater than zero), then there exists a gain a for which the control is 
stable. 

The convergence of the controlled error from the initial state eo to zero is 
governed by Eq. (3.10). It is clear that although small values for a enhance stability 
20 robustness, they also slow down control convergence rate. It is also clear that an accurate 
estimate of the influence coefficient reduces the necessity to be conservative in choosing 
the gain a for stability robustness reasons. Therefore an on-line estimation strategy can 
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be employed to ensure accurate influence coefficient estimates regardless of changing or 
nonlinear dynamics. 

On-Line Estimation of Influence Coefficient 

If the balance correction state is measurable, plant input and output information is 
5 available with which to estimate the influence coefficient after each control iteration. 
The proposed approach for accomplishing the estimation is a recursive exponentially 
weighted averaging method. This method is based on using the current and most recent 
previous measurements to calculate the instantaneous value of the influence coefficient 
using the following formula: 

10 ' c new = Eq.(3.13) 

Previous researchers have utilized non-recursive and recursive least-squares 
estimation using data from every control iteration in past history. This has the effect of 
estimating the average influence coefficient over time without responding well to time- 
varying dynamics. Other work used the initially measured input and output along with 

15 the current measurements compared to the current and most recent previous 
measurements used for this work. In the former estimation technique, if the plant 
dynamics or unbalance disturbance change, or are nonlinear, this method may lead to 
control instability or, at the very least, slower control convergence. Furthermore, no 
provision is made to mitigate the effects of measurement noise. 

20 Assuming sufficient excitation (i.e., sufficient change in the control input wu from 

iteration to iteration), the computed value c new will be non-singular. (A supervisory 
method will be discussed later that "turns off estimation and control in certain situations 
where insufficient excitation persists). However, this estimation formula is also directly 
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subject to the noise present in measuring e and w. Therefore, it may not be desirable to 
throw out all previously estimated information of the influence coefficient. To mitigate 
the effect of measurement noise, an exponentially weighted recursive average can be 
used to update the influence coefficient. The updated influence coefficient used in 
5 computing the next control step is then estimated as 

£ '^-^^i^6 Eq ' (3 - 14) 

Where 0</?<l. The "ft" parameter is a "forgetting factor" that governs the extent of 
control system adaptivity. When the value of "ft" = 0, no weight is given to the newly 
calculated influence coefficient and the control system is the same as the non-adaptive 

10 system discussed previously. When the value of is close to one, the adaptive 
controller places more weight on the most recently estimated influence coefficient and 
tends to ignore previous estimates. 

The forgetting factor /? can be chosen to provide an optimum trade-off between 
estimation sensitivity to measurement noise, and speed of estimation convergence. For 

15 the stationary case after the estimation has substantially converged, the statistical 
variance of the exponentially weighted influence coefficient estimate can be defined as a 
function of the variance of c new . Assuming that the new measurement is independent of 
the past estimate, the variance of the estimate is given by: 

Var[c k ] = y^Var[c nc J. Eq. (3.15) 
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Thus, the estimation variance can be reduced by selecting a low value for fi. For values 
of p close to unity, the estimation variance approaches the variance of the calculated 
parameter c new . This relationship is shown in Figure 3.3a. 

The estimation will "track" the actual influence coefficient c as it varies over 
5 time. The disadvantage of choosing a low forgetting factor f3 is that estimation 
convergence rate is reduced. For the case of a step change in the actual influence 
coefficient to a new constant value, the influence estimation convergence is governed by 
the equation 

E[c k ] = c + (c 0 -c)(l-j3) k Eq.(3.16) 
10 where c is the actual influence coefficient, c 0 is the initial estimate and E[c k ] is the 
expected value of the influence coefficient estimate after k iterations. The number of 
estimation iteration required for 95% convergence of the estimate for varying Rvalues is 
plotted in Figure 3.3b. 

Higher J3 results in faster estimation convergence. For systems that have slowly 
15 varying dynamics and high measurement noise, a low J3 value is best. For highly 
nonlinear or rapidly varying dynamic systems, a larger p value would be necessary for 
fast control convergence. 

Recursive Supervisory Adaptive Control Law 

With the addition of the control gain a and the on-line influence coefficient 
20 estimation, the adaptive single plane recursive control law is given by 

[ e A Eq.(3.17) 
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where c k is given in Eq. (3.14). The block diagram form of the controller is shown in the 
Figure below 

When the magnitude of the change in balance correction vector I w k - Wk-i\ 
approaches zero, the estimation of c k becomes singular. This could lead to "bursts" of 
5 spurious estimates and the balance weights being, commanded to non-ideal positions. 
Though the estimation and control would likely re-converge after a brief "burst" of high 
vibration, such an event is undesirable and could lead temporarily to harmful vibration 
levels. Applicant utilizes a rule-based supervisory strategy. 

The supervisory strategy used for these experiments was to "turn off' the 

10 estimation whenever the change in control vectors - w^i fell below some 
predetermined low level. This level was chosen base on the balance correction resolution 
and measurement noise of the balance correction vector. The control was deactivated 
whenever the vibration error level fell below a preset limit. When deactivated, the 
controller could continue to monitor the vibration error and reactivate the control and 

15 estimation when the error exceeded a certain limit. The limit level included some 
"hysteresis" to prevent annoying on-off cycles of the controller. Thus, the vibration error 
limit at which the controller. Thus, the vibration error limit at which the controller turned 
back on was somewhat higher than the limit at which the controller turned off. A flow 
chart of this supervisory strategy is shown in Figure 3.5. 

20 Such a supervisory control strategy is important for machine tool applications as 

active balance corrections should only be performed prior to, not during, metal cutting. 
The balancing system was configured so that it could receive an "activate" command, and 
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after successful control of vibration error below the limit level, would deactivate itself. 
Vibration error, would be monitored continuously for history logging purposes. 

Single Plane Adaptive Control Experimental Results 

Experimental System Con figuration 
5 The adaptive control law derived above in Eq. (3.17) was implemented for a high- 

speed machining spindle active balancing system. The system was tested on two 
different spindles with similar results. Test spindle 1 was a 23,700 rpm, lkW air-cooled 
spindle. Test spindle 2 was a 1 5,000 rpm, 75k W liquid-cooled spindle. 

A Hall Effect sensor detected a fixed magnetic target on the rotating shaft to 
10 compute spindle rotational speed and to act as a phase reference. The synchronous 
vibration error phasor was computed each control iteration as discussed previously. 
Thirty-two (32) samples of spindle housing radial vibration (measured with a 
piezoelectric accelerometer) were obtained every spindle revolution with sampling 
synchronized to spindle rotation. Data from multiple spindle revolutions (typically six) 
15 were then time-synchronous averaged to provide additional filtering of noise and non- 
synchronous components. 

Balance correction was measured by detecting the passing of a magnet target on 
each balance rotor by stationary Hall Effect sensors. The phase lag of these events 
compared to the shaft phase reference target was calculated to determine the angular 
20 location of each balance rotor. 

Results of Control Nonlinear "Hardening" System 

Results of the single-plane active balancing tests showed the efficacy of the 
adaptive control approach. Test spindle 1 was used to test the adaptive control system 
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robustness in the presence of nonlinear dynamics. The spindle exhibited a "hardening" 
stiffness effect when vibrating at large amplitudes. At the test rotational speed, the 
spindle dynamic stiffness was about two times greater at high unbalance levels than at 
low unbalance levels. For example, on unit of unbalance would cause on unit of 
5 vibration, but two units of unbalance would result in only about 1.5 units of vibration. A 
non-adaptive linear influence coefficient based control algorithm generally may have 
difficulty quickly achieving low vibration in the presence of such a harsh nonlinearity. 
Control stability may also be in question depending on the initial conditions. For the test 
described here, an initial estimate of the influence coefficient magnitude and phase in the 

10 range of previously measured varying values was supplied to the controller. Figure 3.7 
shows how the measured spindle housing vibration at the operating speed 10,000 rpm 
was reduced to an acceptable level by the adaptive controller. 

Figure 3.7 shows that vibration was controlled to below the low limit in two 
balance correction iterations the time required for the vibration data sampling, estimation, 

15 and control computation is evident by noting the time between the end of the first balance 
iteration and the start of the second (observable when the vibration again begins to 
decrease again just before one second). A portion of this time between balance weight 
positioning was allotted as a fixed delay to allow the transient vibration to settle. 

An interesting and surprising side benefit of actively balancing the "hardening" 

20 nonlinear system was the broad frequency band effectiveness of the synchronous 
vibration control. Since the available forcing of the active balancing system can only 
occur synchronous to the machine rotation, active balancing in general can only control 
synchronous vibration. However, for nonlinear systems, large unbalances can excite 
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nonlinear effects that appear at higher harmonic frequencies. Therefore, by controlling 
synchronous vibration, these broadband vibrations can also be significantly diminished. 
This previously unanticipated result is shown in Figure 3.8, a time frequency graph of 
vibration during a single-plane active balancing test. Note that the rotor-synchronous 
5 component is the dominant frequency component of vibration. Other integer and half- 
integer vibration harmonies also appear prior to active balancing. 

After active balancing was completed the synchronous vibration component was 
dramatically reduced. Furthermore, the higher harmonies were also eliminated due to the 
reduction of the synchronous component and the subsequent nonlinear effects. 

10 Tests of Control With Time-Varying Dynamics 

Another important test of the adaptive control system is its robustness to changes 
in plant dynamics. As described previously, rotating machine dynamics can vary over 
time due to changing operating environment, degradation or damage of components, or 
various other reasons. This tendency for rotordynamics to change appreciably highlights 

1 5 the need for a robust adaptive control strategy. 

In the case of high-speed machining spindles, different tool inertias and varying 
machine configurations can result in significant variation of the system dynamics. The 
most typical scenario for time-varying dynamics on a machine tool would be an abrupt 
step change in both machine dynamics and unbalance disturbance after a tool change. 

20 Thus, the "time-varying" characteristic is that the dynamics are unknown and the 
influence coefficient estimate is inaccurate at the beginning of the next active balancing 
cycle. The performance of the single-plane adaptive active balancing system was tested 
for the worst case scenario by providing the control system with an enormous initial 
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estimate of the system influence coefficient. The results of the test with the spindle 
operating about 10,00 rpm are shown in Figure 3.9. 

It is evident that the erroneous initial influence coefficient estimate caused the 
controller to move the balance weights in the "wrong" direction, causing undesired 
5 temporary increase in vibration which could have permanently harmed the tool assembly. 
After each control iteration, however, new input-output data was utilized in a weighted 
average to re-estimate the system unbalance response vector. After four control iterations 
and less than three seconds, the active balancing system controlled vibration to below the 
preset low limit. The controller then had acquired, or "learned", the spindle dynamic 
10 response at that operating speed during the exercising of control. Thus, the active 
balancing system can be used as a probe to measure and track the inherent dynamic 
response of a system. Such data contains useful diagnostic information such as and 
without limitation data associated with the dynamic bearing stiffness and effective 
rotating inertia. 

15 If balancing were again required for the same operating speed and machine 

configuration the control system would be able to reject new disturbances much more 
quickly. Results of a test of the active balancing system after such adaptive "learning" 
are shown in Figure 3.10. 

The test results shown in Figure 3.10 indicate how quickly the single-plane 

20 control system can respond when an accurate estimate of the influence coefficient vector 
is available. 

Results of Control at Multiple Rotational Speeds 



A test of the active balancing system on test spindle 2 resulted in significantly 
decreased vibration at all speeds. The spindle was run at various constant speeds from 
5,000 rpm up to its maximum operating speed of 15,000 rpm and active balancing at each 
speed and is shown in Figure 3.11. 
5 The spindle manufacturer considered 2 mm/sec to be acceptable but less than 1 

mm/sec vibration desirable to protect spindle bearing life. 

During this experiment, the adaptive active balancing system reduced 
synchronous vibration by 79-99%. The control system was able to maintain synchronous 
vibration at or below 0.5 mm/sec even when initial vibration well exceeded acceptable 
10 levels. Background vibration "noise" and the positioning resolution of the balance mass 
actuators determined the minimum achievable controlled vibration level. 

The balance correction states providing minimal synchronous vibration at each 
speed during the tests are shown in Figure 3.12. 

The amount the balance correction magnitude and phase had to change over the speed 
15 range is a measure of how "upstream", or co-located, the vibration control was. If the 
optimal balance correction remained the same over the entire speed range, we would 
know that the balance correction was exciting the same vibrational modes as the 
disturbance unbalance. In this case the control input would be entering the system in the 
same way as the disturbance and the control could be considered to be eliminating the 
20 source of harmful vibration. If the balance correction were required to change 
significantly over the speed range it would be evidence of more "downstream" or non- 
colocated control. That is, the control would be exciting different vibrational modes than 
the disturbance unbalance. The control, therefore, could only eliminate the effect of the 
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disturbance at the sensor, not the disturbance itself. In the case that balance correction 
plane is non-colocated with the disturbance, or the unbalance disturbance is distributed in 
many planes, multiple-plane balancing may be in order. Multiple-plane balancing is 
more capable of "upstream" control than single plane balancing because more vibrational 
5 modes are potentially excited. In the experiments shown above the optimal balance 
correction magnitude and phase only changed slightly above 9,000 rpm. At 7,000 rpm, 
however, balance correction state was significantly different than at other speeds. These 
results indicate that single plane balancing would likely be sufficient at 9,000 rpm and 
above. At 7,000 rpm, results indicate that single plane balancing would likely be 

10 sufficient at 9,000 rpm and above. At 7,000 rpm, additional modes are excited by the 
residual unbalance disturbance which may require multiple plane balancing to control. 

A significantly useful side benefit of the active balancing system is that it can be 
used as a probe to measure and track the unbalance response of a system. Such system 
identification information contains useful data about the dynamic bearing stiffness and 

15 effective rotating mass. For the experiment discussed above, the unbalance sensitivity 
was stored at each speed every time the active balancing system was enabled to control 
vibration. This data was obtained for two spindle system configurations to allow for 
comparisons. A long tool and short tool were separately inserted into the spindle and 
active balancing control enabled for successive rotational speeds. Figure 3.13 and Figure 

20 3.14 below show the measured unbalance sensitivity amplitude and phase angle for the 
two tool configurations. 

These data were obtained by activating the automatic balancing system at each 
speed and storing the last influence coefficient estimate used during active balancing. 
- • • 
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Thus the unbalance response shown also indicates the control authority magnitude of the 
active balancing system on spindle vibration. 

These last tests were conducted with the spindle mounted on a massive test bench. 
Since the tool mass was considerable less than the overall vibrating system mass, the 
5 unbalance response did not change dramatically when the different tools were inserted. 
One can observe in Figure 3.12 that, as one might expect, the critical speed close to 
11,000 rpm did shift to a slightly higher frequency with the smaller tool inserted. The 
peak response also decreased slightly. 

Stability of Adaptive Influence Coefficient Control 

10 The single plane adaptive influence coefficient recursive control law was earlier 

derived. The controller 10 commands the balance correction at control iteration £+1 to 
the state given in Eq. (3.17) by the use of "a" and "P" control parameters. The manner in 
which this is accomplished is further shown below. 

For the case in which the actual plant influence coefficient and the disturbance 

15 remain constant during control convergence, the overall adaptively controlled error 
response can be derived. Such a scenario is representative of the typical machine tool 
active balancing application. After each tool change, balancing is performed at the new 
operating speed at which cutting is to be performed. The relevant dynamics can, in 
general, be expected to vary significantly from tool change to tool change because of 

20 different spindle operating speeds, toolholder inertias, and geometric configurations of 
the machine structure. However once the spindle is accelerated to the new operating 
speed, the unbalance disturbance and machine dynamics do not typically change during 
the few seconds in which active balancing is performed. 
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The adaptively controlled vibration error response can be formulated using the 
control and estimation equations previously derived. This formulation is based on the 
assumptions that the actual influence coefficient c and unbalance disturbance d do not 
change during control convergence, and that "sufficient excitation" exists to allow non- 
singular estimation at each control iteration. When the control has converged well enough 
that the balance correction is no longer changing, the estimation must be "turned off. 
After such convergence the active balancing control response will be the same as non- 
adaptive control method. To analyze the control response during convergence, we can 
substitute the expected value of the influence coefficient estimate from Eq. (3.16) into the 
control of law Eq. (3.17). The adaptively controlled vibration error response is then given 
by the recursive relationship 

f \ 

ac ■ 

e *« = 1_ 777^ „vi R\ k Y k Eq.(4.1) 

V c + (c 0 - c)(l - P) ) 

The adaptively controlled vibration error can then be described by the infinite product 



n 



1- aC Eq. (4.2) 



c + (c 0 -c)(l-/^ 

Recall that the estimation forgetting factors defined such that 0 < fi < 1. 
Therefore, regardless of the initial influence estimate c 0 , the argument of the infinite 
product of Eq. (4.2) converges (again assuming sufficient excitation) in the limit such 
that: 



( ac 
lim 1 



= l-a ' Eq.(4.3) 



c + (c 0 -c){\-/J) J ) 

Because the control gain a is defined such that 0<a <1, there will exist a control 
iteration p (p < oo) such that when j = p the magnitude of the argument will be less than 

23 



one. Although the infinite product argument magnitude will eventually fall below one 
and thereafter cause the product to converge, the value of may continue to increase 
until the p th iteration. The corresponding vibration error e p will reach the value 

v*fl('- e+ft -Tx.-J *™ 

5 After the p th iteration, the error magnitude will begin to decrease. Thus the magnitude of 
the error at the p ili iteration depends on the magnitude of the estimation error (c 0 - c) and 

the value of /?. An example of this is shown in Figure 4.1 . 

Note that the error magnitude in Figure 4.1 continues to increase until the infinite 
product argument falls below one at the third control iteration. In this case the iteration p 
10 =3 and the maximum worst case error magnitude is computed from Eq. (4.4) to be \e p \ 
= 7.31. 

Because it occurs before the infinite product argument converges to 1-oc, the value 
p will always be finite. Therefore, the maximum controlled error must be bounded. For 
any arbitrary finite value e p and any complex term a such that \a\ < 1, the infinite 
15 product shown in Eq. (4.5) will exponentially converge to zero as k -> oo.' 

k 

=^n« Bq. (4.5) 

- J=P 

The significance of this result is that, assuming sufficient excitation and that the plant 
influence coefficient and disturbance remain constant during control convergence, the 
infinite product of Eq. (4.2), and hence the adaptive control law of Eq. (3.17) is stable 
20 regardless of the initial influence coefficient estimate. 
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The stability analysis above is based on the assumption that neither the unbalance 
disturbance nor the machine dynamics changed during control convergence. However, it 
is reasonable to extend the stability conclusion to certain cases of time-varying 
disturbances and influence coefficients. Simplistically speaking, after initial 
5 convergence, the on-line estimation would "track" changes in the machine dynamics with 
exponential convergence as long as the influence coefficient varied less than 2/3 times the 
previous influence coefficient each control iteration. This is because the estimation will 
respond a fraction J3 times the distance to the actual value during each iteration. 
Similarly, if the unbalance disturbance changes less than 2 a times the previous value 

10 during each iteration, the control would still exhibit exponential convergence because the 
control would get closer to the optimal value each iteration. These comments are 
somewhat simplistic, however, because they neglect the more complicated issue of the 
combined effect on estimation of disturbances and influence coefficients that change 
from iteration to iteration. However, the statements do provide a general idea of what 

15 rates of change the disturbance and dynamics will still allow stable control. 

The only caveat in the stability condition proven above is that there is no way of 
telling how high the error signal gets before it begins to converge to zero. Choosing 
of=j3=1 will cause faster adaptation and control convergence. However, the higher the 
value of a, the worse the temporary vibration error will become during control 

20 convergence. Thus, although the control system is stable an erroneous initial influence 
coefficient estimate will always present a trade-off between speed of control 
convergence, and the amount of temporarily high vibration "overshoot". 
Parameter Auto-Tuning Adaptive Control 
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If one were able to vary the adaptive parameters a and ft from control iteration to 
control iteration, one might be able to eliminate the trade-off between control 
convergence rate and the worst case temporary vibration. Such a parameter "auto- 
tuning" method would, therefore, allow the best of both worlds during control 

5 convergence. When the influence estimate was accurate (i.e. the infinite product 
argument in Eq. (4.2) was close to the value of 1-a) a value of a close to one would 
provide the most rapid control convergence. When the influence estimate was not 
accurate, a low value of a would provide more cautious control and thereby limit the 
worst case vibration error "overshoot". A method of automatically selecting the 

10 parameters a and ft according to this concept is hereafter described. 

Since the actual influence coefficient during control is unknown, it is, 
unfortunately, impossible to know how accurate the influence estimate is. It is, however, 
possible to measure the convergence of the influence estimate. One could accomplish 
this by for estimation scheme discussed here by comparing the difference between the 

15 instantaneously measured influence coefficient, and the exponentially weighted averaged 
estimate. A parameter can be defined as a normalized measure of influence coefficient 
estimation convergence error at iteration k such that 



where c new and c k _ x are defined in Eq. (3.13) and Eq. (3.14) respectively. The symbol 
20 "*" represents the complex-conjugate. Assuming that the actual influence coefficient and 
unbalance disturbance did not change, this estimation convergence error measure will be 
unbiased. The parameter <% will eventually converge to zero as the influence estimate 




Eq. (4.6) 
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converges according to Eq. (3.16) and will be calculated by controller 22 for each control 
iteration. 

To implement the automatic adjustment of parameters a and {3 during each 
control iteration, the following functions are proposed. 

5 a k =\-e { ' ]/M) Eq.(4.7) 

J3 k = e^ 1 "*^ Eq. (4.8) 

where a k and fi k are the control gain and estimation forgetting factors respectively to be 
used at each control iteration k by the controller 22 and where 7j a and rjp are arbitrary 
real scaling factors to allow flexibility in shaping the response of each parameter. Eq. 

10 (4.7) and Eq. (4.8) ensure that when the influence estimation error is high, the control 
gain ak is low. Examples of acceptable scaling factors include the value of 1 for each 
factor. This will tend to limit and/or substantially reduce the worst case transitory 
vibration error. As shown in Figure 4.2 when the estimation error is low, the control gain 
pCk approaches one, ensuring fast control convergence. Conversely, the estimation 

15 weighting parameter fi k is close to unity when estimation error is high, ensuring fast 
estimation response, is close to unity when estimation error is high, ensuring fast 
estimation response, is low when the estimation is more accurate, reducing the effect 
of measurement noise. 

Note that if estimation error was ever exactly zero, the exponents in Eq. (4.7) and 

20 Eq. (4.8) would both become singular. The exponential, however, would approach zero 
in this case. For the practical implementation of the method, a very small positive 
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number can be added to the denominator in the exponents during computation to ensure 
numerical robustness. 

Using the function for j8 k given in Eq. (4.8), the auto-tuning estimation equation is 
given by 



where a* is given by Eq. (4.7) and c k is computed from Eq. (4.9) and where the applied 

balance relates to and/or comprises the movement of the rotors 28, 30 to a 

10 position defined by the value of "w£+f (e.g. a position providing the level or amount of 
unbalance correction represented by the value of " Wk+ 
Automatic Supervisory Limit Selection 

A supervisory control strategy is necessary for many machine tool applications 
because active balance corrections should only be performed prior to, not during, metal 

15 cutting. The balancing system should be configured so that it could receive an "activate" 
command and, after successful balancing, would deactivate itself and signal the machine 
controller to begin machining. If the balancing system were not deactivated, the 
nonstationary synchronous vibration signal likely to be present during machining could 
lead to the active balancing device constantly "hunting" to find the optimal balance 

20 correction. Balance adjustments based on the vibration during cutting could potentially 



5 




Eq. (4.9) 



The auto-tuning adaptive control law is subsequently 




Eq.(4.10) 
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cause an increase in unbalance and harmful vibration. Furthermore, such continuous 
operation requires more power to the balance actuators and could lead to undo wear and 
tear on the device. To alleviate these problems, a supervisory strategy is in order to 
prevent excessive balance adjustments and estimation bursting. 
5 A practical consideration for end users of such an active balancing system is how 

to set the vibration error limits to activate and deactivate control appropriately. One 
would want to set them as low as possible to ensure that low vibration levels were 
achieved. However, measurement noise, variation in machine dynamics and resolution 
limitations of the active balancing device constrain just how low the vibration error can 

10 be controlled. Because of these issues, some specialized engineering knowledge must be 
applied in selecting supervisory limit vales. Such specialized knowledge is not always 
available to the typical end-user of the active balancing system. Furthermore, when 
machine speed or setup changes, the limit levels would require re-adjustment. Therefore, 
an automated method of selecting control limits is proposed here that incorporates the 

15 considerations mentioned above. 

For active balancing devices with discrete states such as the one used for this 
research, the worst case correction resolution can be defined. This resolution, in 
combination with the estimate of the system influence coefficient, can be used to define 
the low vibration error limit at which control will "deactivate". This limit can be defined 

20 as 

^'■= Kk««" M Eq. (4.11) 

where c k is the estimated influence coefficient a the current control iteration k w re s is the 
worst case balance correction resolution (in the same unbalance units used in the 
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influence coefficient), 7] e is a unit-less scaling factor (simply for added flexibility in 
shaping response) and s k is the influence coefficient estimation error defined in Eq. (4.6). 
The exponential term is included to ensure that the control does not prematurely 
deactivate because of an erroneous influence coefficient estimate. When the influence 
5 coefficient estimation is converged (^ is small) the low control limit will be the vibration 
error expected at the worst case resolution. 

Using single activation-deactivation vibration limit may result in spurious control 
activation if random noise causes the measured vibration error to jump back above the 
limit. To prevent this, an additional "high" limit can be defined to incorporate a 

10 hysteresis band into the supervisory control. Thus once the control is deactivated when 
the vibration falls below the low limit, it is not reactivated unless the vibration exceeds a 
higher limit as shown in Figure 4.3. 

For the automated limit selection method used here, after controlling vibration 
error to below the "low" limit of Eq. (4.1 1), control will not be reactivated until vibration 

15 error magnitude exceeds a "high" limit defined as 

e h ighk =3a e +e Wi . Eq. (4.12) 

where cr e is the standard deviation of the vibration error magnitude measurement. This 
standard, deviation could be continuously measured whenever the control system was 
deactivated and simply monitoring the idling spindle vibration error. 

20 The low limit of Eq. (4.1 1) was chosen so that, in the worst case, the mean value 

of vibration error magnitude will be exactly at the low limit. In the most cases, the mean 
error magnitude will be less than the low limit. The measured error magnitude will 
fluctuate about the mean value because of measurement noise. Assuming that the 
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measurement noise is Gaussian, even in the worst case when the error magnitude mean 
were stationary at the low limit, there would be a 99.7% probability that the measured 
error magnitude would not exceed the high limit of Eq. (4.12). therefore, by 
automatically defining the high limit this way, there is a high probability that control 

5 would not be spuriously reactivated purely because of measurement noise. 

By using Eq. (4.11) and Eq. (4.12) to automatically select supervisory control 
limits, no specialized user knowledge or user input is required. Furthermore, these 
automatic selection criteria take into account specific plant dynamics, active balance 
correction resolution and vibration error measurement noise. 

10 Experimental Results 

An active balancing device was mounted on a conventional and commercially 
available high-speed grinding spindle and the auto-tuning adaptive control tested for 
various cases of unknown step changes in the unbalance disturbance and influence 
coefficient. 

15 Particularly, the test spindle was a Fischer Model No. MFW- 1240 grease- 

lubricated liquid cooled 10 kW high- speed grinding spindle. This would be a typical 
state-of-the-art spindle for CBN grinding applications. The spindle allowed for safe 
testing at the relatively high rotational speed of 20, 000 rpm. 

Auto-Tuning Control Performance Compared With Typical Fixed- Parameter 
20 Control With Erroneous Influence Coefficient Estimate 

Tests were performed to compare the auto-tuning control with the conventional 

fixed-parameter adaptive control developed. Figure 4.5 shows the comparative results 

for the condition of an erroneous initial estimate of the influence coefficient. The initial 

25 influence estimate used was approximately 180° out of phase with the actual influence 
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coefficient. The conventional fixed adaptive parameters were set to typical values of a = 
0.9and/?=.08. 

Note that the erroneous initial influence coefficient estimate causes both 
controllers to temporarily increase the vibration error magnitude. However, the auto- 
5 tuning controller allows a much lower worst error magnitude while converging as fast, or 
even faster than the conventional adaptive controller. The corresponding filtered 
synchronous vibration error values and adaptive parameter values are shown in Figure 
4.6. 

Comparison With "Cautious* Fixed-Parameter Control With Erroneous 
10 Influence Coefficient Estimate 

To limit the worst case error magnitude to that of the auto-tuning controller, the 

conventional controller must use an initial control gain of a =0.1. Since the parameter 

remains constant during conventional control convergence, however, significantly 

15 reduced speed of control convergence results. Figure 4.7 shows the results of such a 
"cautious" conventional adaptive control approach compared again to the auto-tuning 
controller performance. 

The corresponding filtered synchronous vibration error values are shown in 
Figure 4.8. Setting the control gain to a low value allowed the conventional adaptive 

20 control to limit the worst case overshoot vibration error. However, the convergence rate 

was significantly slowed. 

Comparative Performance With Typical Fixed-Parameter Control With Accurate 

Influence Coefficient Estimate 

25 When the initial estimate is accurate, conventional adaptive influence coefficient 

control performance is typically very good. The worst case error is very low and control 



32 



convergence is rapid. Figure 4.9 shows the comparative performance of the auto-tuning 
and conventional adaptive control when an accurate influence estimate is available. 

Figure 4.10 shows the corresponding filtered synchronous vibration error values 
and adaptive parameters during the control convergence. 
5 In the case of an accurate influence estimate, the auto-tuning control performs just 

as well as the conventional adaptive control. In fact the only situation in which the auto- 
tuning control will not perform at least as well as the conventional control is when the 
conventional control can converge in one step. Since a low initial value of control gain 
ctk is used in the auto-tuning control, the control will never converge in only one step. It 
10 would often be able to converge in two steps in the same conditions, however. 

Summary of Experimental Auto-Tuned Adaptive Control Performance 
The adaptive parameter auto-tuning functions introduced above enhance 
conventional adaptive influence control by limiting the worst case temporary vibration 
15 error while allowing the same (or better) control convergence rate. Furthermore, because 
parameters are varied automatically during each control iteration, no user setup is 
required. This is typically a distinct advantage because no special operator training is 
required to setup and use the vibration control system for each individual machine or 
environmental condition. 
20 Automatic Supervisory Limit Calculation Experimental Results 

The automatic supervisory limit selection method was enabled during the 
experimental testing. Figure 4.11 shows the automatically calculated supervisory limit 
settings from the experiment shown in Figure 4.6. 
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Control was deactivated once vibration error magnitude was controlled below the 
low limit at about 1.0 seconds elapsed time. Control was never reactivated because 
vibration error magnitude never exceeded the high limit. The standard deviation of the 
vibration error magnitude measurement noise during this test was 0.0029 g's. The active 
5 balancing device used in the test had two stepper-motor type balance rotors each with 60 
detent increments per revolution. The worst case balance correction resolution for this 
configuration is nln^ums (i. e., 5.2%) of the maximum balance correction capacity of the 
device. 

In summary, the previous description shows that Applicant's invention 
10 substantially eliminates performance trade-offs inherent in the fixed-parameter adaptive 
control and automatically incorporates engineering knowledge to simplify the active 
balancing system operation for the end-user. 

The adaptive influence coefficient method allows active balancing control without 
a priori modeling or experimentation to determine plant dynamics. A stability analysis 
15 was presented for the fixed-parameter adaptive control in the typical case where the 
synchronous disturbance and plant influence coefficient control was shown to be stable in 
these cases regardless of the initial influence coefficient estimate. However, the worst 
case temporary vibration during control adaptation could become quite large for 
erroneous influence coefficient estimates. A trade-off exists between speed of control 
20 convergence and the magnitude of this worst case error. 

An automatic tuning method was therefore presented that allows the adaptive 
control parameters to be adjusted each control iteration to limit the low worst case error 
magnitude while still providing fast control convergence. The parameters were adjusted 
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each control iteration based on an instantaneous measure of the influence coefficient 
estimation accuracy. When estimation was inaccurate, the control become less 
aggressive and the estimation placed more weight on the instantaneously measured 
influence coefficient. When estimation error was low, the adaptive parameters were 
5 adjusted to provide more aggressive control and to place more weight on the long-term 
averaged estimation. 

Supervisory control is necessary to "turn off control so that machining 
operations can proceed once vibration is controlled below an acceptable limit. 
Significant specialized engineering knowledge is typically required to set supervisory 

10 vibration error limits for each application. Automatic supervisory limit selection criteria 
were defined to eliminate the need for such specialized end user input. The vibration 
error limits for enabling and disabling control were defined based on functions of 
estimate of machine dynamics, estimation accuracy, active balance mass actuator 
resolution and vibration measurement noise. 

15 Experimental results showed that the "auto-tuned" controller provided control 

response speed comparable to, or better than, the conventional adaptive control while 
maintaining significantly lower the worst case vibration error magnitudes. The 
experiments also illustrated the effectiveness of the auto-tuning supervisory control with 
integrated automatic error limit calculation. The multiple plane dithering aspect of 

20 Applicant's invention will now be described. 

Derivation of Optimal Control Law 

The multiple-plane influence matrix-based active balancing control law was 
previously discussed. For square influence matrices [C], the control law is again 
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presented here with the estimated influence matrix [c] used instead of the actual 
influence matrix, which is not necessarily known. 
Wk+i = Wk- [c]' l {E}i 

For the case of more error sensors than balance planes, the influence matrix is not 
5 square. Because the influence matrix does not exist in this case, the control law of Eq. 
(5.1) cannot be used. Active balancing, therefore, cannot in general drive all the error 
signals to zero. Researchers involved in off-line balancing have noted that least squares 
solution (or "pseudo-inverse") of influence coefficient balancing can be implemented for 
non-square influence matrices. Such a control scheme allows the minimization of the 
10 sum of the squares of error signal residuals. In the event that not every error sensor 
output were deemed as important as every other, a weighted least squares control law can 
be implemented. A real positive semi-definite n x n (where n is the number of error 
sensors) diagonal penalty matrix [Q] can be specified to give the desired relative 
weighting of the error sensors: 



15 [Q] = 



q . 0 

12 



Eq. (5.2) 



0 q„_ 

where q, is the relative weight of the i' h error sensor. Similar real positive semi-definite m 
x m (where m is the number of active balance planes) weighting matrices [R] and [S\ can 
be defined to penalize control effort, and control rate of change respectively. An 
objective function J({E}u+u {W}k+\, {W}k, [Q], [R], [5]) can then be defined such that 
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J=\{e)I x [q]{e} m +\{w)Ur]{w} m +\({w} m -\w} k )\s]{{w} k+} -{w} k ) 

Eq.(5.3) 

Where the "*" symbol denotes the complex -conjugate transpose operator. The optimal 
control problem then consists of commanding the balance weight vector for the next 
5 control iteration { W}k+j so as to minimize the objective function J. 

Assuming that the disturbance vector {D} does not change over one control 
iteration, the error vector relationship of Eq. (2.5) can be applied for two control 
iterations and the results combined to yield the relationship 

{E} k+} = [C]({W}m ~{ W}k) + {E) h Eq. (5.4) 

10 Substituting Eq. (5.4) into Eq. (5.3) allows evaluation of the objective function J 

using only the updated control vector { W} k +\ and measured values {E} k and { W) k \ 

, , . Eq. (5.5) 

To find the stationary points of the objective function J with respect to the 
updated control vector {W} k +u we can take the corresponding partial derivative, set it 

15 equal to zero and solve for the control vector. This stationary point is guaranteed to be a 
minimum if the objective function J is truly a quadratic function of {W} k+ \. To ensure 
that J is a quadratic function, at least one of the matrices [Q] 9 [R] 9 and [S] must be 
positive definite. Furthermore, if only [Q] is positive -definite, then the matrix [C] must 
also be full rank (rank m in this case). To ensure that the influence matrix [C] is full 

20 rank, care must be taken to install the active balancing devices in appropriate planes. 
These are: 1) avoid placing them at nodal locations of any vibrational modes to be 
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controlled; and 2) make sure that the influences from each balance plane on the error 
vector are independent for all speeds at which the active system will be operated. 

Using matrix calculus techniques, the equation for minimal objective function is 
found to be 

Eq. (5.6) 

The optimal control vector update { W}k+\ can be solved for algebraically as follows: 
([C\'[Q][C] + [R] + [S]){W}m = aQ\Q][C\ + [S\){W}k- [Q\Q]{E} k ) 
Eq.(5.7) 

io = + + Eq. 

(5.8) • 

At steady-state, assuming that control has converged, the optimal control vector is then 
given by 

{m«,o pt =-([Q\Q][Q+[R]y\cf[Q]{D} Eq.(5.9) 
1 5 and the subsequent steady- state error given by 

{E}* opt = ({i m „} - [Qm\Q][QHm\Q\Q]){D} Eq, (5.10) 

where [I nxn ] is the n x n identity matrix. Note again that the practical realization of the 

optimal control law in Eq. (5.8) must use the estimated influence matrix [c] instead of 

the generally unknown or varying actual influence matrix. 
20 Performance and Stability Analysis of Non- Adaptive Optimal Control 

The [R] matrix allows penalizing the control effort at each active balance plane, 
which could be desirable in certain cases. For example, if the source of residual 
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unbalance were not co-located with the active balance planes, the optimal balance 
correction could be different for different operational speeds. The optimal balance 
correction at operating speed may sometimes cause harmful vibration at another slower 
speed if an emergency shutdown resulted in deceleration through a critical speed faster 
5 than the active balancing control could track. By conservatively limiting the control 
effort using a nonzero [R] matrix, higher vibration levels at operating speed could be 
traded off for a reduced possibility of causing harmful vibration during an emergency 
shutdown of the rotating machine. 

The [5] matrix provides for penalizing the speed of control response. This can 

10 have the benefit of potentially allowing enough time for operator intervention in the case 
of any sort of malfunction of the active balancing system. 

Both the [R] and [S] matrices provide opportunities for cautious control that is 
often required in very conservative industries utilizing turbomachinery in critical 
continuous processes. Furthermore, both matrices can have the added benefit of 

15 enhancing the stability margin, of non-adaptive control. By limiting the control inputs or 
slowing down the control response, an implementation that might originally be unstable 
due to a bad influence matrix estimate can be made stable. An analysis of the optimal 
control stability follows. 

Replacing the actual influence matrix in the optimal control law of Eq. (5.8) with 

20 the estimated influence matrix and substituting in the error vector response relationship of 
Eq. (2.5), we obtain the recursive relationship 

Eq.(5.11) 
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which leads to 
5 Eq.(5.12) 

For stable control, the control input {W}k+\ must converge to a constant value. Only the 
first term multiplied by {W} k is germane to the stability question. The second term 
multiplying {D} is constant and only affects the control vector to which the steady-state 
control converges. It then follows that the optimal active balancing control is stable for a 

10. constant disturbance vector {D} and constant influence matrices [C] and [c] if and only 

if 



([C]^ 



<1 Eq.(5.13) 



where the symbol a signifies the maximum singular value. The effect of the matrix [R] 

on stability is evident in Eq. (5.13). Regardless of the estimation error ([c] - [C\\ a 

15 large enough [R] matrix will render the control "stable" by effectively preventing the 
control effort from increasing unbounded. The effect of the control rate penalty matrix [S\ 
on stability is perhaps not quite so clear from Eq. (5.13). The stability criterion can be 
expanded again as 

|c] - [el<?] + w +[«])"■ ([£]' [e][c] - [c] - [ejc] + [s]) 

20 Note that since [c]*[£>][c] is positive semi-definite, if only the real part of 

[c] *[£7][C] is positive definite, and if both [c] and [C] are full rank (rank m), then there 



< 1 Eq. (5.14) 
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will exist a positive semi-definite real matrix [S] for which the stability criterion of Eq. 
(5.13) is met. This is analogous to the situation described for the single-plane control of 
Chapter 3 where the influence coefficient estimate is at least within ±90° of the actual 
value. The single-plane control, in that case, could be stabilized using a smaller control 
5 gain a. 

Assuming stable control, and assuming that the disturance vector {D} and actual 
influence coefficient [C] are constant during convergence, it can be shown that the 
balance correction vector will converge at steady state to 

=-([c]\qIc] + [R]Y[c}'[q]{D} Eq. (5.15) 
10 The corresponding steady-state error vector will be 

'■{*}. =([/«H^ Eq.(5.16) 
It is apparent from Eq. (5.15) that even if control is stable, if the influence matrix 
estimate [c] does not equal the actual influence matrix [C], then the steady-state balance 

correction vector will not necessarily equal the optimal correction vector given in Eq. 

15 (5.16) will not equal the error vector of Eq. (5.10) and the control will not be the same as 
with a perfectly accurate influence matrix estimate. This is one key limitation of the non- 
adaptive control approach. 

If one applies non-zero[7?] or [S] penalty matrices purely to enhance stability 
robustness of the non-adaptive control, certain performance trade-offs become necessary. 

20 By inspecting Eq. (5.15), it is evident that as [R] approaches infinity, the steady-state 
balance correction approaches zero and the steady-state controlled error vector 
approaches the disturbance {£>}. Thus the stability robustness associated with a large [R] 
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matrix comes at the expense of reduced steady-state control weighted least-squares 
control performance. If the estimated influence matrix element phase angles fall within 
90° of the actual element phase angles such that Re([c ]*[£?] [C]) is not positive definitely 
then the [S] matrix can be increased to enhance stability margin. This comes at the 
5 expense of slowing down the control convergence rate. The steady-state results in this 
case also could also could be non-optimal in the same sense as mentioned above. 

The stability margin trade-offs and non-optimal performance associated with non- 
adaptive control using erroneous influence matrix estimates can all be remedied by 
ensuring an accurate influence matrix estimate through on-line system identification. The 
10 next section details the implementation of such a strategy for multiple-plane active 
balancing systems. 

On-Line Estimation of Influence Matrix 

To enhance stability robustness and improve performance optimality of the 
influence matrix-based control, the adaptive control developed in Chapter 3 can be 
1.5 extended for multiple-plane applications. The multiple-plane estimation also utilizes the 
current and most recent previous sets of measurements to calculate the instantaneous 
' value of the influence matrix. This instantaneous value is then recursively averaged with 
the previous estimate. 

Upon initial startup, or after an unknown step change in the disturbance vector 
20 {D} and influence matrix [C], we again assume these vectors do not change during 
control convergence. Eq. (2.5) for m successive control iterations (where m = the number 
of balance planes) can then be combined to calculate the influence matrix by the 
following formula: 
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[q^[A£][A W\ A Eq.(5.17) 
where [A£] and [AW] are matrices containing the changes in correction vectors and error 
vectors respectively after m successive control iterations such that 



Eq. (5.18) 
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Eq. (5.19) 



Note that the [AJV\ matrix must be at least (m x m) and non-singular to compute [C] wew . 
Thus measured data corresponding to (m +1) independent balance correction states must 

10 be obtained to re-calculate the influence coefficient matrix. 

Assuming sufficient excitation (i.e., sufficiently independent changes in, the 
control vector {W) k form iteration to iteration), the computed value [C] mw will be 
nonsingular. (A supervisory method will be discussed later that "turns off" estimation 
and control in certain situations where insufficient excitation persists). As in the single- 

15 plane estimation method detailed earlier, the effect of measurement noise can be 
mitigated using an exponentially weighted recursive average to update the influence 
matrix estimate. The updated influence matrix used in computing the next control step is 
estimated as 

[^L - fr-^FL +4A£L(AFT]- , A Eq. (5.20) 
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where again the "forgetting factor" is defined as 0 ;<JS < 1 and can be used to govern the 
extent of control system adaptivity. 

Multiple Plane Dithering Method 

As discussed previously, to re-compute the influence matrix each iteration, m 
5 independent balance correction movements are required to form a non singular "Delta" 
matrix as shown in Eq. (5.18). However, as the on-line estimation of the influence matrix 
converges to the actual value, the controller would tend to output nonindependent 
correction vector commands at each successive control iteration. To alleviate this 
problem, the control must be "dithered" somehow to provide m independent correction 

10 vectors. The simplest way to accomplish this would be to move the correction weights 
only one balance plane at a time (while recording the error signal vectors after each 
balance weight movement). This approach, however may cause error signals at some 
sensor locations to become temporarily worse after moving the balance weights at only 
one correction plane. The potential for worsening vibration depends on the magnitude of 

1 5 cross-influence present. Cross-influence is when off-diagonal elements of the influence 
matrix are close to the same magnitude as diagonal elements. The more cross-influence 
in the system, the more chance for temporarily generating worsened vibration. To ensure 
that all error sensor output amplitudes decrease, all correction weights could be moved 
simultaneously (assuming an accurate influence matrix estimate). However, this would 

20 not provide independent measurements with which to calculate an updated influence 
matrix estimate. By creating a matrix of independent basis vectors one could vary the 
relative independence of each balance sub-iteration and thus have some control over the 
"singularity" of [AW]. Consider the symmetric matrix [7] such that 
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Eq. (5,21) 



where m is the number of balance planes, and y (0 < y < 1) is an output "dithering" 
parameter that defines the independence of the balance correction vectors for each control 
5 sub-iteration. Note that the sum of all the elements in each row of the basis matrix [7] is 
unity. Thus at the end of all the sub-iterations the correction weight in a given plane will 
be fully at the commanded {{W}k+\ state. The commanded balance weight vectors for 
the p th sub-iteration of the control iteration (k+l) would be 

{W} k+hp = {W} k+i ^ 1+ ([l mxm ]{T} p )({W} k+] ~{W}^ Eq.(5.22) 

10 for integer pe [0, m], and {W}k+\,o initializes as {W}k- Umxm] indicates the m x m identity 
matrix and {T} p indicates the p th column of the basis matrix [7]. Note that for y=l the 
matrix [7] is exactly the identity matrix. Subsequently each sub-iteration would involve 
moving one balance plane correction at a time for maximum correction vector 
independence and maximum influence matrix estimation signal-to-noise ratio. At the 

15 other extreme for y= 0, each sub-iteration would involve moving each balance correction 
plane simultaneously the same amount. This would minimize the worst temporary error 
signal amplitude (assuming an accurate influence matrix). This would also, however, 
result in a singular [AW] matrix and render estimation of the influence matrix impossible. 
In practice, the y parameter can be specified to provide the desired trade-off between 
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multiple plane influence coefficient estimation signal-to-noise ratio and worst-case error 
signal amplitude during active balancing. 

Recursive Optimal Adaptive Control Law 

A recursive adaptive implementation of the optimal control law of Eq. (5.8) can 
5 be obtained by the addition of the on-line estimation given in Eq. (5.20) and the 
"dithering" method governed by the basis matrix of Eq. (5.21). The recursive optimal 
adaptive control law is then given by the equation 

Eq. (5.23) 

10 

where [cj* is the estimated influence matrix computed using Eq. (5.20). Each control 

iteration k is divided into m sub-iterations over which the correction vector is "dithered" 
according to Eq. (5.2 1 ) to ensure non-singular estimation. 

If the correction vector {W} k does not change from iteration to iteration, the 

15 influence matrix estimation will also become singular. After estimation "bursting", the 
estimation and control would be expected to converge again. However, to prevent such 
"bursting" and the associated high controlled error levels, a supervisory strategy can be 
implemented to "turn off estimation when the correction vectors do not change 
significantly (i.e. [A W\k becomes singular) and to turn off control when the error vector is 

20 controlled below a pre-determined level. A flow chart for the supervisory 
implementation of the optimal adaptive control is shown in Figure 5.1. 
Experimental Results of Optimal Multiple-Plane Control 
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The optimal adaptive control law was implemented for a two-plane active 
balancing system on a laboratory flexible rotor test rig. The test rig consisted of a 0.65" 
(16,5mm) shaft supported over a span of approximately 30 inches (762 mm) on two ball 
bearings and driven through a narrow "quill-shaft" coupled to a direct current type motor. 
5 Active balancing devices were mounted to the shaft at approximately the third-span 
locations. Each balance actuator and disk assembly weighed approximately 10 lbs. (4.5 
kg). Eddy current proximity probes mounted close to, and outboard of the two disk 
locations were used to measure radial shaft deflection in two orthogonal (i. e., x and y) 
directions. 

10 The first two critical speeds of the rotating test rig were measured to be 

approximately 1 lOOrpm and 4360rpm. The predicted mode shapes for these two critical 
speeds are shown in Figure 5.3. 

To test the active balancing system for a plant with fairly "rich" rotordynamics 
the rotor was run at about 3,100 rpm during the experiments. Since this speed was in 

15 between the first two critical speeds, both mode shapes contributed significantly to 
dynamic response. Figure 5.4 shows the predicted forced-response deflected shape of the 
rotor due to a 0.34 oz-in unbalance at the drive-end (left end) active balancing plane. 
This unbalance represented the maximum correction capacity of the active balancing 
device. Figure 5.4 could be used to predict the influence coefficient at 3,100rpm from the 

20 drive-end active balancing device to the two sensors. The combination of contributions 
of the two mode shapes at 3,100 rpm provides for a relatively interesting response 
situation. Note that unbalance at the drive end balance has a much greater influence on 
shaft response at the outboard sensor than on the drive-end sensor. Since the rotor 
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configuration is highly symmetrical, the outboard balance plane also has the greatest 
influence on the drive-end sensor. The actual experimental results matched this 
prediction fairly well in a qualitative sense. The average measured influence matrix 
between both balance planes and four sensors at 3,100 rpm is presented in Table 5.1. 
5 Single-Plane Versus Multiple-Plane Active Balancing Results 

Because two modes contribute significantly to response at 3,100 rpm, one active 
balance plane would not in general be sufficient to balance the rotor. Figure 5.5 shows 
the results of active balancing using one balance plane (drive-end) compared with results 
using two balance planes. The adaptive control law was used in both cases and all 

10 optimal penalty weighting matrices had the same values [Q] = [7], and [R] = [S] = [0]. 

Note that the initial vibrations for the one and two-plane tests were not the same 
because the initial active balance correction in the single-plane balancing experiment was 
not zero (residual unbalance for both tests was the same). It is evident from Figure 5.5 
that the single-plane optimal control does not match the performance of the dual-plane 

15 control. In fact, the single-plane control is only able to improve a small amount over the 
initial vibration error. This is an indication that: 1) residual unbalance was exciting more 
than one mode; and 2) the disturbance vector lays mostly outside the span of the single- 
plane "controllability" space. The latter means that the residual unbalance was not 
concentrated at the active balance plane and did not excite the two modes in the same 

20 way as the active balance correction. 

Weighted Output Optimal Adaptive Control 

As noted above, there may be occasions where balance correction should be 
restricted to allow for safe traversing of different vibration modes in the event of an 
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emergency shutdown of the machinery. In this case the control effort penalty matrix [R] 
could be increased to penalize correction magnitude. Figure 5.6 compares adaptive 
control results for both and non-zero [R] matrices. 

Note that the penalized steady-state control effort (balance correction magnitude) 
5 is reduced to about one half of the non-penalized control effort. This optimal trade-off 
obviously also results in higher steady-state vibration error magnitudes. 

Partially Output Dithered Adaptive Optimal Control 

The control output dithering required to ensure non-singular estimation could 
potentially cause temporarily high vibration error values even when a good influence 

10 estimate is available. This effect will be worst when there is a high level of "cross- 
coupling" between balance planes. That is, the off-diagonal terms of the influence matrix 
are significant relative to the diagonal terms. Figure 5.7 compares the results of adaptive 
control shown in Figure 5.6 using full output dithering (y=l) versus the performance 
using reduced output dithering (y=0.25). Although the influence matrix in this case was 

1 5 not necessarily "highly coupled", one can still observe a difference in the first control 
iteration. 

Results of Fully Output Dithered Adaptive Optimal Control of Lightly Cross- 
Coupled System 

20 Note that the vibration in the first iteration of fully dithered output actually 

increased over the initial vibration. By updating only one active balancing plane at a 
time, balance correction does not move in a "straight line" toward the optimal state' The 
excursions of the correction state about the optimal state can potentially cause such 
worsening vibration before the other balance plane is updated. Using only partial output 
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dithering allows both planes to be partially updated in the same control iteration, thereby 
reducing the "excursions" of the balance correction state from the optimal state. The 
lower dithering parameter number causes the balance correction vectors to proceed more 
"smoothly" toward the optimal states. This is illustrated by the continuously decreasing 
5 objective function value for the partial output dithered control in the figure above. 

Comparison Of Non-Adaptive And Adaptive Optimal Control With Somewhat 
Inaccurate Influence Matrix Estimate 
A significant test of the optimal control methods was to provide an erroneous 
influence matrix estimate to the controller. Figure 5.8 shows a comparison of optimal 
10 control results four control scenarios when the initial influence matrix estimate is 
erroneous. The influence matrix estimate for the test was approximately 0.6 times the 
magnitude and rotated 70° from the actual matrix shown in Table 5.1. The first plot of 
Figure 5.8(a) shows the unstable performance resulting from traditional non-adaptive 
control under this condition. Figure 5.8(b) shows how the non-adaptive control was 
15 stabilized using a non-zero [R] matrix. Thirdly, Figure 5.8(c) illustrates non-adaptive 
control stabilization using a non-zero [S] matrix. Figure 5.8(d) shows the results of 
adaptive control. 

In each case the error sensors all received equal weighting and the [R] and [S] 
matrices were either zero matrices if not noted or were diagonal with the value shown. 
20 Because of the inaccurate influence matrix estimate, the conventional non-adaptive 
weighted least-squares control is unstable as shown in Figure 5.8(a). The active 
balancing devices were both saturated (outputting maximum correction possible), 
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providing the only limit to the vibration increasing further. Such vibration values 
represent damaging and potentially dangerous levels. 
Dual-Plane Optimal Active Balancing Control Results for Somewhat Inaccurate 

Initial Influence Matrix Estimate 

5 Adding a non-zero penalty on control effort to the objective performance function 

has the positive effect of producing stable control as seen in Figure 5.8(b). The only 
downside is that the steady-state vibration levels are not especially low. The use of a 
non-zero control rate of change penalty term in the objective performance function also 
results in stable control as illustrated in Figure 5.8(c). It is evident, however, that the 

10 control convergence rate is not ideal. Figure 5.8(d) shows that adaptive control not only 
provides stable control, but converges relatively quickly. This is in spite of the fact that 
the output must be dithered to ensure non-singular estimation. The "overshoot" of the. 
vibration levels for the; adaptive control case was caused by the erroneous influence 
matrix initial estimate. After one control iteration, however, the influence matrix 

15 estimate was corrected and the vibration subsequently was reduced rapidly to low levels. 

The corresponding objective function values for the results of Figure 5.8 are given 
in Figure 5.9. Also shown in Figure 5.9 are estimates of the stability criterion of Eq. 
(5.13) for each optimal control scenario. This criterion could not, in general be measured 
directly as the actual influence matrix is unknown. For Figure 5.9, the average influence 

20 matrix measured during experimentation (shown in Table 5.1) was used as the "actual" 
influence matrix [C] in Eq. (5.13). ' 

For stable control, the stability criterion must be below the value one. Note that 
the unstable traditional non-adaptive control has a stability criterion value of close to two. 
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By adding nonzero penalty matrices [R] and [5], the stability criterion value drops below 
one and the control converges. The stability criterion for the adaptive control begins at 
the same unstable value as the traditional non-adaptive control. This corresponds to the 
temporarily worse vibration in the first control iteration. However, as the on-line 
5 estimation converges to the actual influence matrix, the adaptive stability criterion value 
falls below one and eventually falls close to zero. The result is that the control converges 
rapidly. 

The non-adaptive control is rendered stable using the non-zero [5] matrix. 
Furthermore, the performance in terms of weighted least-squares error is better than the 

10 case of using a non-zero [R] matrix. However, it is important to note that the steady-state 
performance using the non-zero [S] matrix was still sub-optimal because of the inaccurate 
influence matrix estimate. The minimum steady-state objective function value attained 
with the non-adaptive control was 0.16. Though much lower than the uncontrolled value, 
this was almost 2.5 times greater than the minimum objective function value of 0.0665 

15 attained using the adaptive control. Thus it is clear that adaptive control not only 

provides fast and stable control convergence, but also provides more optimal steady-state 

performance as predicted by observing Eq.(5. 1 5) and Eq. (5.1 6). 

Comparative Optimal Control Results for Highly Inaccurate Influence Matrix 

Estimate 

20 When the influence coefficient estimate has any elements more than 90° rotated 

from the corresponding element actual influence matrix, a non-zero [5] matrix cannot 
stabilize the non-adaptive control. By choosing an appropriate non-zero [R] matrix, the 
unstable outputs of the non-adaptive control can be effectively bounded. Such control 
will not have optimal results, however. In fact, the vibration cannot be reduced but will 
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only be prevented form increasing out of control. Figure 5.10 shows the results of using 
an influence matrix estimate whose elements are all rotated approximately 135° from the 
actual matrix elements. 

The non adaptive control experiments in Figure 5.10(a) was concluded with only 
5 a few control iterations because the balance corrections were diverging and had saturated 
at full balance correction capacity (correcting in the "wrong" direction). By increasing 
the control effort penalty, although the vibration becomes worse after control, the steady- 
state balance correction capacity does not saturate at full capacity. Hence the steady-state 
vibration seen in Figure 5.10(b) is not as bad as the traditional non-adaptive case. The 

10 test shown in Figure 5.10(b) was concluded after the active balance corrections had 
achieved steady state and were no longer changing. The control results shown in Figure 
5.10(c) again illustrate the effectiveness of the adaptive control strategy even when the 
influence matrix is initially estimated poorly enough to cause unstable non-adaptive 
control. After initially increasing the vibration, the adaptively controlled vibration 

15 quickly converges to very low levels. 

It is to be understood that the inventions are not limited to the. exact construction 
or methods illustrated and described above. Rather, various changes may be made 
without departing from the spirit and scope of the inventions as defined and as are set 
forth in the following claims. 

20 
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